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INTRODUCTION 



X 



Stars form in cluste r s and associ ations (e.g. ICarpenterll2000l ; 
lLada fc Ladal 120031 ; lLadal l20lrf >. and their spatial distri- 
bution appears to be smooth and contin uous, i.e. there 
is no evidence of bi-modal star formation (|Bressert et al.l 
|2010| ). However, the morp hologies of individual clusters ap - 
pear to vary significantly (|Cartwright fc Whitworthll2004 ). 
with some regions displaying a high amount of substruc- 
ture (e.g. Taurus and Chamaeleon) and others displaying 
a smooth, centrally concentrated morphology (e.g. p Ophi- 
uchus and IC348). This suggests that either star formation 
in clusters is not universal and may depend on the local envi- 
ronment, or that all clusters form with the same morphology 
which is later altered by dynamical interactions. 

There are several methods to look for structure in 
cluster s, and even clus t ers t hemselves. The study by 



ABSTRACT 

We determine the distribution of stellar surface densities, £, from models of static and 
dynamically evolving star clusters with different morphologies, including both radially 
smooth and substructured clusters. We find that the £ distribution is degenerate, in 
the sense that many different cluster morphologies (smooth or substructured) produce 
similar cumulative distributions. However, when used in tandem with a measure of 
structure, such as the Q-parameter, the current spatial and dynamical state of a star 
cluster can be inferred. The effect of cluster dynamics on the £ distribution and the 
Q-parameter is investigated using ./V-body simulations and we find that, depending on 
the assumed initial conditions, the £ distribution can rapidly evolve from high to low 
densities in less than 5 Myr. This suggests that the £ distribution can only be used 
to assess the current density of a star forming region, and provides little information 
on its initial density. However, if the £ distribution is used together with the Q- 
parameter, then information on the amount of substructure can be used as a proxy to 
infer the amount of dynamical evolution that has taken place. Substructure is erased 
quickly through dynamics, which can disrupt binary star systems and planets, as well 
as facilitate dynamical mass segregation. Therefore, dynamical processing in young 
star forming regions could still be significant, even without currently observed high 
densities. 

Key words: stars: formation - kinematics and dynamics - open clusters and associ- 
ations: general - methods: numerical 



and can be use d to search for clusters , or "clustering " in 
crowded fields (|Gutermuth et alj 120091 ; ISchmeial [iuTir i. in 
tandem with the mean separation between stars, to define 
the Q-parameter. The Q-parameter quantifies whether the 
cluster is substructured or radially concentrated, and to 
what extent, for a given morphology. 



ICartwright fe Whitworthl (|2004T ) used a minimum spanning 
tree analysis, which has also been u sed to quantify the 
amount of mass segregation in clusters IjAllison et al. 2009a) 



The stellar surface density distribution has al so been 
used t o study the structure of star forming regions. iLarsonl 
i| 19951 ) showed that the structure in Taurus was hierarchi- 
cal, but broke down at the b inary regime ( corresponding to 
the local Jeans length), and ISimonl |l997l ) found the same 
result for Ophiuchus and the Trape zium cluster. However, 
iBate. Clarke fc McCaughreanl (|l998l ) repeated this analysis 
for the Trapezium cluster and found that the break between 
random stars and binaries does not co rrespond to the Jeans 
length in this cluster. Furthermore, iKraus fc Hillenbrand! 
(2008) found that the break between random stars and bi- 
naries was llOOOau in Upper Sco, compared to 17000au 
in Taurus, and suggested the break is a function of the age 
(and initial density) of the cluster, and is caused by dynam- 
ical evolution. 
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Recently, the distribution of stellar surface densities (£) 
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was also used t o evaluate the d efiniti on of what constitutes 
a star cluster. iBressert et al. (2010) examined a volume- 
limited sample of nearby star- forming regions and calculated 
the distribution of E in each region. IBressert et all (|2010h 
then compared the cumulative distribution of E for all the 
regions (apart from the central region of the Orion Nebular 
Cluster (ONC), where the determination of E was not pos- 
sible for every star due to crowding) to various definitions of 
clustering, based on stellar surface density thresholds. Due 
to a lack of obvious bi-modality in the distribution, they 
concluded that stars in the local solar neighbourhood form 
in a continuous spatial distribution. 

In this paper, we determine the effectiveness of the 
surface density distribution, E, in tandem with the Q- 
parameter, for providing information on the dynamical state 
of star clusters. We begin by analysing static star clusters 
with various morphologies to look for variations in the shape 
of the E distribution. We then dynamically evolve radially 
smooth Plummer sphere clusters, and substructured frac- 
tal clusters, and determine the E distribution and the Q- 
parameter as a function of time. We describe the set-up of 
the clusters in Section [5] we present our results in Section [3] 
we provide a discussion in Section [4] and we summarize and 
conclude in Section 



2 METHOD 

In this section we describe the measurement of stellar surface 
density, E, and the Q-parameter, which is used to quantify 
the amount of substructure in star cluster. We then describe 
the Monte Carlo set-up of static star clusters and associa- 
tions before describing the set-up of iV-body simulations in 
which we measure E and the Q-parameter in star clusters 
as a function of time. 



2.1 The stellar surface density distribution, E 

The recent work bv lBressert et al.l l|2010l ) measured the stel- 
lar surface density, E, of companions in a volume-limited 
sample of star forming regions. They adopt the following 
formula l|Casertano fc Hut 1 19851 ): 



N ■ 



(1) 



where N is the N th nearest neighbour and Djv is the pro- 
jected distance to that nearest neighbour. The value of E 
can be sensitive to the choice of N; for example a low 
N value would cause an artificially high E value if the 
binary fraction were high. A high N value would fail to 
detect localised pockets of high density, which can have 
importan t implications for the evolution o f substructured 
cluste rs (|Parker. Goodwin fc Allison! l201ll ) . IBressert et al.l 
(|2010h adopt N = 7, which is also the value we choose. 



2.2 The Q- parameter 

The Q-parameter wa s pioneered by 

ICartwright fc Whitworthl |2004l and combines the nor- 
malised mean edge length of the minimum spanning tree of 



all the stars in the cluster, m, with the normalised correla- 
tion length between all stars in the cluster, s. The level of 
substructure is determined by the following equation: 

a-s. m 

A substructured cluster has Q < 0.8, whereas a smooth, cen- 
trally concentrated cluster has Q > 0.8. The Q-parameter 
has the advantage of being independent of the density of 
the star forming region, and purely measures the level 
of substructure present. The original formulation of the 
Q-parameter assumes the cluster is spherical, but can 
be altered to take into account the effects of elongation 
|Bastian et al.ll2009l ). 



2.3 Cluster morphologies 

We set up clusters with four different morphologies. Two 
are radially smoot h, centrally con centrated morpholo gies - 
a Plu mmer sphere (jPlummeri 191 ll ) , and a King profile ijKind 
1962). Both are used extensively as the initial conditions of 
JV-body simulations and assume a relaxed, spherical sym- 
metry. We would therefore expect that these clusters would 
display a smooth, continuous distribution in stellar surface 
density. 

The remaining two morphologies are intended to mim- 
ick the observations of young, dynamically unevolved 
star forming r egions (like the majority of those in the 
IBressert et al.l (|2010l) sample), which show a high level of 
subst r ucture (e . g. iGutermuth et al . 2009; S anchez fc Alfarol 
120091 : ISchmeial and references therein). Firstly, we 

adopt the fractal distribution. Secondly we create 'associ- 
ations', with several 'nodes' of star formation. 

For each morphology, we create a cluster with N = 1000 
stars, and then find the physical size scale (i.e. cluster ra- 
dius) for that morphology w hich has the best ove rlap with 
the observed distribution in IBressert et"al] l|2010l ). To this 
end, we simply choose the cluster radius so that the median 
E best corresponds to the median value in the Bressert et 
al. distribution (Eyso = 22 stars pc -2 ). The random num- 
ber seed used to generate the positions of stars in our model 
clusters results in some variation around this median, but 
we are examining the morphologies of these clusters, rather 
than attempting to produce an exact fit to the observed 
data. Note also that our chosen sample of stars is smaller 
than the 3857 used to make the cumulative distribution in 
IBressert et~ai1 l|2010l ). The E distribution in Bressert et al. is 
the sum of 12 different star forming regions, with a median 
membership of N ~ 100, but is dominated by the Orion 
complex (2696 objects). All we aim to investigate here is 
which, if any, single morphology could feasibly reproduce 
the observed E distribution. 



2.3.1 Plummer spheres 

Plummer spheres l|Plummerlll91ll ) were first used to describe 
the radial profiles of Globular clusters, and are used exten- 
sively to model the initial mass distribut ion in star clusters 
in iV-body simulations (e.g. lKroupall2008l ). due to their sim- 
plicity. The mass enclosed in radius r is given by the follow- 
ing equation: 
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where M c i us is the total mass of the cluster, and r p i is 
the characteristic 'Plummer radius'. Formally, the Plummer 
sphere extends to an infinite distance, but the radius within 
which half the mass resides, fi/2, is given by ri/2 = 1.305r p i 
(|Kroupall200i ). 



a sphere from the cube and are then randomly removed (so 
maintaining the fractal dimension) until the required num- 
ber of children is left. These children then become stars in 
the star cluster. 

The fractal is therefore described by two numbers, the 
fractal dimension D, and the radius r>. It is not clear 
whether star form ation actually gives a fractal d istribution 
in nature (see e.g. lElmeereen fc Elmegreenll200l1 ). but it is 
the most convenient method for creating substructure. 



2.3.2 King profiles 

King profiles (|Kind 1 19621 ) are a more realistic description 
of of the radial profile of Globular clusters. The King pro- 
file is used extensively in A-body simulations, but also to 
fit the profiles of observed clus ters, for example the ONC 
( Hille nbrand fc Hartmaimll 19981 ). A core radius, r c (and cen- 
tral surface density, fo) are defined in the model, and the 
point at which the radius of the cluster is tidally truncated, 
rt. The surface density, / at a distance r from the centre of 
the cluster is given by the following relation: 



fa 



1 + (r/r c ) 2 



(4) 



King profiles can be described by a central 'concentration 
parameter' , Wo (with W o oc logio (r/rt)), and the truncation 
radius, rt (King 1966). If a King profile has a large value 
(> 3) of Wo then it has an extended envelope, whereas a 
King profile with a low Wo is said to just have a central 
core. We will detail our choice of Wo and rt in the following 
Section. 



2.3.3 Fractals 

The fractal distribution has the advantage that the substruc- 
ture is described by one parameter, the fractal dimension, 
D. In three dimensions, a cluster with fractal dimension 
D — 1.6 is very clumpy, whereas a cluster with D = 3.0 
is a uniform sphere. We set our fractal clusters up ac- 
cordin g to the method outlined in iGoodwin fc Whitworthl 
(2004). This begins by defining a cube of side iVdiv (we 
adopt iVdiv = 2.0 throughout), inside which the fractal is 
built. A first-generation parent is placed at the centre of 
the cube, which then spawns A| iv subcubes, each contain- 
ing a first-generation child at its centre. The fractal is then 
built by determining which of the children themselves be- 
come parents, and spawn their own offspring. This is deter- 
mined by the fractal dimension, D, where the probability 



that the child becomes a parent is given by 7V ( 



(D-3) 



For a 



lower fractal dimension, fewer children mature and the final 
distribution contains more substructure. Any children that 
do not become parents in a given step are removed, along 
with all of their parents. A small amount of noise is then 
added to the positions of the remaining children, preventing 
the cluster from having a gridded appearance, and the chil- 
dren then become parents of the next generation. Each new 
parent then spawns A| iv second-generation children in iV| iv 
sub-subcubes, with each second-generation child having a 



(D-3) 



probability of becoming a second-generation parent. 



This process is repeated until there are substantially more 
children than required. The children are pruned to produce 



2.3.4 Associations 

We also create a substructured association, by randomly 
placing 5 'nodes' of star formation in three dimensions. From 
each node, we distribute stars in a spherical coordinate sys- 
tem, defining a characteristic length scale, or 'compactness', 
C which we divide the distance from the centre of each node 
iVby: 



N x + ^-^sin#cos(? 

N y H — — sin0sm<p 
O 



N : 



XR „ 
+ —^-cosO, 

O 



(5) 
(6) 
(7) 



where A is a random number between and 1, 9 £ [0, 7r], 
4> G [0, 2ir] and R is the radius of the association. 

For each cluster morphology we determine E for each 
star and make a cumulative distribution of the E values. We 
then determine the Q-parameter for each cluster. Note that 
we 'observe' our clusters in the x-y plane; for each cluster 
we have checked that the results do not depend on viewing 
angle (any difference is of order a few per cent - a similar 
effect to changing the random number seed used to generate 
the positions). We detail our choices of cluster parameters 
(e.g. ri/2, D, Wo, etc.) in Section I3TT1 



2.4 Dynamical evolution 

We study effects of dynamical evolution on the distribution 
of E and the Q-parameter for two types of cluster morphol- 
ogy - a radially smooth and centrally concentrated Plum- 
mer sphere, and a substructured fractal. Both the Plummer 
sphere and the fractal comprise of N = 1 000 singl e stars 
with masses drawn from a 2-part power law lKroupal (|2002l ) 
IMF with the following parameters: 



N(M) 



M~ 2,3 



m Q < M/Mq ^ m 3 
mi < M/Mq ^ m 2 



(8) 



where mo = 0.1 Mq, mi = 0.5 Mq, and m 2 = 50 Mq. Our 
choice of m 2 is due to our adoption of an ONC-like cluster 
to simulate; the most massive star in the ONC, 6 Ori C has 
a system mass of 45-50 M Q (|Kraus et al.ll2007l . 120091 ). We 
do not include brown dwarfs in the simulations. We do not 
include primordial binaries, because they do not affect the 
determination of E. 
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Clusters are unlikely to remain gravitationally bound 
if they undergo an early phase of gas expulsiorrj 
llTutukovll 19781; lLada et al.lll984f ; [Goodwin fc Bastianll200rJ 
iBastian et al.ll2008h . and we model this effect by setting half 
of our simulations to be supervirial, which causes them to 
immediately begin to expand. We also set up clusters in 
virial equilibrium (static), and subvirial (collapsing). 



2.4-1 Plummer sphere clusters 

The stars have positions and velocities chosen accord- 
ing to the prescription from a Plummer sphere in 
lAarseth. Henon fc Wielenl (|l974T ). We set up clusters in two 
different 'states'; one has a half-mass radius, ri/ 2 = 0.8 pc 
and is in virial equilibrium at the start of the simulations 
(a virial ratio of Q v ir = 0.5, where Qvir = T/\Q\ and T and 
Q. are the total kinetic and potential energy of the stars, re- 
spectively). The second cluster has a much smaller half-mass 
radius {r\/2 = 0.1 pc) and is initially supervirial (i.e. initially 
expanding), with Q v i r = 1.5). 



Table 1. A summary of the different cluster properties adopted 
for the TV-body simulations. The values in the columns are: the 
number of stars in each cluster (Altars), the morphology of the 
cluster (either a Plummer sphere or fractal) , the initial virial ratio 
of the cluster (Q v ir)j the initial radius of the fractal, (rp), or the 
initial half-mass radius of the Plummer sphere, (r^ m). 



iVstars 


Morphology 


Qvir 


r F or n/ 2 


1000 


Plummer sphere 


0.5 


0.8 pc 


1000 


Plummer sphere 


1.5 


0.1 pc 


1000 


Fractal 


0.3 


1 pc 


1000 


Fractal 


1.5 


1 pc 



3 RESULTS 

In this section we show the E distributions for various clus- 
ter morphologies, and calculate the Q-parameter (in 2D), 
before showing the effects of dynamical evolution in radially 
smooth, and substructured clusters on both the E distribu- 
tion and the Q-parameter. 



•-4. 



Fractal clusters 



We set up the fracta ls with an initial radius of lpc accord- 
ing to the method in lGoodwin fe Whitworthl (|2004l V As de- 
tailed in Section 12.31 the level of substructure in a fractal 
is determined by just one number, the fractal dimension, 
D. We set up highly substructured fractals, with D — 1.6. 
As with the Plummer sphere clusters, we draw 1000 stars 
from the IMF in Equation [8] and distribute them randomly 
in the fractal. The velocities of the 'parents' in the fractal 
(see Section r2.3.3[) are drawn from a Gaussian of mean zero, 
and the 'children' inherit their parent's velocity plus a ran- 
dom component that decreases with each generation of the 
fractal. This means that nearby stars have similar velocities 
to their neighbours, but distant stars can have very different 
velocitiefl 

Finally, we scale the velocities of the stars in the fractal 
to the desired virial ratio; in one set of simulations we 
choose a subvirial (cool) virial ratio, Q v i r — 0.3, and in the 
second set of simulations we again adopt a supervirial ratio 
(Qvir = 1.5). 

We run 10 suites of each simulation, identical apart 
from the random number seed used to initialise the posi- 
tions, masses and velocities of the stars. The simulations are 
evolved for 1 Myr using the kira integrator in the Starlab 
package (e.g. iPortegies Zwart et alJll999l . l200lh . We do not 
include stellar evolution in the simulations. Details of the 
four simulations are given in Table [1] 



1 We note that the recent study bv lKruiissen et"aT. I <2012h finds 
that gas expulsion may not have a strong effect on the evolution 
of clusters in hydrodynamical simulations; see also the discussion 
of recent ob servational evid ence on the relative importance of gas 
expulsion in Bastian ( 2 0111 '). 

2 A discussion of the potential effects of non-correlated velocities 
on the dynamical evolu t ion o f the fractal clusters is provided in 
iGoodwin fc Whitwortt] \2004) . 



3.1 Static clusters 

For each of the four morphologies, we perform a two-sample 
Kolmogorov-Smirnov (KS) test between the E distribution 
for the model cluste r and the observed distribution from 
iBressert et all (|2010f i. If the test returns a p-value of less 
than 0.01 (arbitrarily chosen to be "low") we reject the null 
hypothesis that there is no difference between the two dis- 
tributions. A summary of the KS p-values for all four mor- 
phologies is given in Table(2] along with the level of structure 
measured by the Q-parameter. 



3.1.1 Plummer spheres 

A Plummer sphere with half-mass radius 2.0 pc is shown in 
Fig. 1(a) with the corresponding E distribution shown in 
Fig. 1(b) At first glance, the Plummer sphere appears to 
provide a reasonable fit to the observed E distribution in 



local star forming regions, and in particular it reproduces 
the tail of the distribution well. 

However, a KS test on the Plummer sph ere data and 
the observational data in IBressert et all (|201Ch returns a p- 
value of < 10 -8 and we reject the null hypothesis that the 
two datasets are drawn from the same distribution. 

The the level of structure measured by the Q- 
parameter for this Plummer sphere is Q = 1.41, con- 
sistent with a smooth, cen trally concentrated profile 
|Cartwright fc Whitworthll2004h . 

3.1.2 King profiles 

A King profile with a radius of 10 pc and a concentration pa- 
rameter Wo = 12 (w hich corresponds to an extended outer 
envelope, lKindll966l ) is shown in Fig. This shows a better 
fit to the observational data than the Plummer sphere, with 
the only deviation occurring at high surface densities. 

A KS test returns a p-value of 0.02 and for this mor- 
phology we cannot reject the null hypothesis that the two 
datasets are drawn from the same distribution. 
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Figure 1. A Plummer sphere with half-mass radius rj/2 = 2.0 pc, and the corresponding cumulative di stribution of st e llar s urface 
densities, S. The S distribution is shown by the solid line, and for comparison we show the data from iBressert et al. I ll2010l l with 
the dashed line. The value of was chosen such that the median £ corresponds to the median value in the observational data 

(Syso = 22 stars pc -2 ). 




Figure 2. A King profile cluster, with a radius, rt of 10 pc and a concentration parameter, Wo = 12, and the corresponding cumulative 
distribution of stella r surface densities, S. The S distribution is shown by the solid line, and for comparison we show the data from 
IBressert et al. with the dashed line. The values of rt and Wo were chosen such that the median S corresponds to the median 

value in the observational data (Eyso = 22 stars pc -2 ). 



We find that Q = 1.14, indicating that the King profile 
is not as centrally concentrated as the Plummer sphere. This 
is due to our adoption of the King concentration parameter 
of Wo = 12, which produces an extended outer envelope 
around the cluster core. 



3.1.3 Fractals 



In Fig. 3(a) we show a fractal distribution (with fractal di- 



mension D — 1.6) and the corresponding E distribution in 
Fig. 3(b) The fractal has an overall scale of 20 pc. The me- 



dian v alue of S for th i s frac tal agrees with the median value 
in the IBressert et al.l (|20T0h data, but the distribution of S 
is much narrower. 



The KS p-value is < 10 and we reject the null hy- 
pothesis that these datasets could be drawn from the same 
distribution. 



The Q-parameter for this fractal is Q — 0.20, indicating 
a high level of substructure. 
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Figure 3. A fractal distribution, with fractal dimension D = 1.6 and an overall radius of 20 pc, and the corresponding cumulative 
distribution of stella r surface densities, E. The E distribution is shown by the solid line, and for comparison we show the data from 
iBressert et al, with the dashed line. The radius was chosen such that the median E corresponds to the median value in the 

observational data (Eyso = 22 stars pc -2 ). 




-20 20 0.01 0.1 1 10 100 1000 

x/pc I [stars/pc 2 ] 

(a) (b) 

Figure 4. An association, comprised of five 'nodes' with a cluster radius of 25 pc, and the corresponding cumulative distribution of 
stellar surface densities, E. The E distribution is shown by the solid line, and for comparison we show the data from IBressert et all 
with the dashed line. The radius was chosen such that the median E corresponds to the median value in the observational data 
(Eyso = 22 stars pc -2 ). 



3.1.4 Associations 

In Fig. [4] we show an association with radius R — 25 pc, 
with five nodes placed randomly in three dimensions. The 
compactness of each node, C, is 7.5. As with the fractal, 
the median value of E agrees well with the mean value from 
the observational sample; however, and like the fractal, the 
distribution is not wide enough. 

The KS p-value is < 10~~ 5 and we reject the null hy- 
pothesis that these datasets could be drawn from the same 
distribution. 



The calculation of the Q-parameter for this association 
gives Q = 0.32, indicating it is substructured, but not to the 
same extent as the fractal. 

3.1.5 Possible biases 

In principle there are several biases that could affect the 
determination of E in star forming regions. Firstly, the E 
distribution could be affected by the erroneous inclusion of 
foreground or background stars. Depending on the density 
of the star forming region compared to the background, the 
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Figure 5. The effect of background stars on the distribution of S. In panel (a) we show a cluster with size 10 pc and a nodal 'compactness' 
C = 7.5, against a uniform background distribution of 1000 stars. The histogram of S is shown in panel (b) and the corresponding 
cum ulative distribution of E is shown in panel (c). The £ distribution is shown by the solid line, and for comparison we show the data 
from Brcsscrt ct al. (2010) with the dashed line. In panel (d) we reduce the compactness of the association nodes, so that C = 2.0, and 
we retain the spatial density of the background stars. The histogram of £ and the corresponding cumulative distribution are shown in 
panels (e) and (f), respectively. 



Table 2. A summary of the four different morphologies of our 
static clusters. We list the 2D Q-parameter for each morphology, 
and the KS test p- value between the S distribution of e ach cluster 
and the observational data from Brcsscrt et al. | <201Cft . 



Morphology 


KS p-value 


Q-parameter 


Plummer sphere 


3 X 10~ 9 


1.41 


King Profile 


2 X 10~ 2 


1.14 


Fractal 


8 x 10- 36 


0.20 


Association 


4 x 10~ 6 


0.32 



inclusion of field stars shows an obvious signature. In Fig. [5] 
we show the same association set-up in Fig. [4] with the same 
compactness, C = 7.5 but with a smaller cluster size (lOpc). 
We then randomly place a further 1000 stars within the 
field of view and determine the distribution of E including 
all stars. The histogram of E clearly shows bimodality in 
the distribution (Fig. |5(b)[ ) and in Fig. |5(c)| we show the 
cumulative distribution, which displays a prominent kink in 
the transition from the cluster nodes to the background. 

We now change the C parameter of the nodes in the 
association to C = 2.0. This has the effect of dispersing 
the stars in the nodes so that they are less distinct against 



the background (see Fig. 5(d) |. The bimodality is no longer 
apparent in the histogram (Fig. 5(e) l (although it is still 



distinguishable by eye) and the cumulative distribution of 



E reverts to a smooth line (Fig. 5(f) I. 

Secondly, the determination of E could be affected by 
"missing" stars, hidden by dust extinction. We examine the 
effects of extinction on a flux limited sample for a Plummer 
sphere cluster by applying a simple extinction power law to 
the cluster. We reduce the level of extinction from A v = 100 
at the cluster centre to A v = at a distance of 2.0r 1/ / 2 , 
where rwa is the half-mass radius, using the following for- 
mula: 



A v (r) = 100 



1.5 n 



/2 



(9) 



where r is the position of the star with respect to the cluster 
centre and \r\ is the modulus of its vector. To account for 
projection effects along the line of sight, we double A v (r) 
if the z-component of the vector r is negative. This en- 
compasses the densest region of the Plummer sphere (see 
Fig. 6(c) I. We have varied the power- law exponent between 



3 and 7 and find very little difference to the results. 

If a star has a mass m < 0.2 Mq and has an extinc- 
tion A v (r) > 25, then it is 'hidden', and not counted in the 
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Figure 6. Plummer spheres with half-mass radii r\ji = l.Opc, and the respective £ cumulative distributions before removal of extincted 
stars [(a) and (b)], and after [(c) and (d)]. The extent of the dust clou d (2,0 r 1 /?) is shown by the circle in panel (c). The £ distribution 
is shown by the solid line, and for comparison we show the data from lBressert et alj l|201fj ) with the dashed line. 



determination of E. We find that 25 per cent of stars are 
hidden if we apply the above constraints. 

The results of extinction on the distribution of E are 
shown in Fig. [6] In panel (a) we show the original Plummer 
sphere (with half-mass radius rj/g = l.Opc), before the ex- 
tincted stars are removed. In panel (b) we sho w the distribu- 
tion o f E by the solid line, with the data from lBressert et al.l 
<|2010h shown for comparison. In panel (c) we show the clus- 
ter with the extincted stars removed (and the extent of the 
dust cloud shown by the circle) and the corresponding E 
distribution in panel (d). The effect of extinction is to move 
the distribution to the left, indicating that a star forming re- 
gion could appear less dense due to 'hidden stars'. However, 
we note that the level of extinction required to produce this 
shift is drastic, and in practice is unlikely to bias observa- 
tional studies. 

Similarly, the determination of the Q-parameter is not 
strongly affected by the removal of extincted stars. For a 



Plummer sphere we found Q = 1.41, which is reduced 
to Q = 1.36 after the removal of 25 per cent of the 
stars in Fig. [6(a)| T his is consistent with the findings of 
iBastian et al.nf2009l ). who find that the Q-parameter de- 
creases by 0.07 ± 0.03 from the true value when stars in a 
sample are hidden by extinction. 



3.2 Dynamical evolution 

In this section we show the cumulative distribution of E 
for clusters that undergo dynamical evolution. We calculate 
the cumulative distribution of E before the cluster evolves 
(i.e. at Myr), then at 1 Myr, 5 Myr and 10 Myr. We first ex- 
amine the evolution of smooth, spherical Plummer spheres 
in two dynamical states; virial equilibrium (Q v ir = 0.5) and 
supervirial, or 'warm' (Q v ir = 1.5). We then apply the same 
analysis to a highly substructured fractal, both subvirial 
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I [stars/pc 2 ] I [stars/pc 2 ] 

(a) Q v i r = 0.5 (virial equilibrium) (b) Q v i T = 1.5 (expanding) 

Figure 7. The evolution of X for Plummer spheres. In (a), we show an N = 1000 Plummer sphere initially in virial equilibrium (Q v i r with 
a half-mass radius of 0.8 pc. In (b) we show the evolution of a 'Hot' Plummer sphere, half-mass radius = O.lpc, virial ratio Qvir = 1.5. 
The data from Bressert et al are shown in red, and the distribution of X for the clusters at OMyr, lMyr, 5Myr and 10 Myr are shown 
by the solid, dashed, dash-dotted and dotted lines, respectively. 



(undergoing cool collapse with Q v i r = 0.3) and supervirial 
'warm' (Qvir = 1.5). 

Additionally, we calculate the Q-parameter and the evo- 
lution of the virial ratio as a function of time for the clusters 
that are initially substructured. When analysing dynamical 
simulations, the Q-parameter can be artificially high due to 
extreme outliers (R.J. Allison, priv. comm.), so we remove 
the outer 2 per cent of stars in the iV-body snapshot from 
its determination. 



3.2.1 Plummer spheres 



In this second scenario, the E distribution for the clus- 
ter before dynamical evolution lies far to the right of the 



In Fig. |7(a)| we show the cumulative distribution of E for 
a Plummer sphere cluster in virial equilibrium, where the 
virial ratio Q v i r = 0.5. The half-mass radius of this cluster 
is initially ri/2 = 0.8 pc, and the E distributions at OMyr, 
1 Myr, 5 Myr and 10 Myr are shown by the solid, dashed, 
dot-dashed and dott ed lines, respectively. The data from 
I Bressert et al.1 l|2010f ) are shown by the solid red line in this 
plot. Although the cluster expands over its lifetime, and the 
E distribution moves to lower values, it does not change in 
the first 1 Myr. 

Alternatively, star clusters could undergo a period 
of rapid expansion due to (almost) instantaneous gas 
expulsion by the first supernovae in the clust er (e.g. 
iBastian fe Goodwinll200rj ; lGoodwin fc Bastianll2006l . though 



Kruijssen et al 



2012) 



For example, the velocity dis- 
persion of th e Orion Nebula C luster is estimated to 
be 4.3 kms" 1 (|Qlczak et al.1 120081 ), although iFuresz et al] 
(2008) find a lower value of 3.1 kms" 1 , which is higher than 
if the cluster were in virial equil ibrium and suggests that 
the c luster is expanding (see also lGieles fc Portegies ZwartJ 
l201ll ). We model this by setting the virial ratio to be 'su- 
pervirial' at the start of the simulation. We show the evo- 
lution of an initially very dense (7^/2 = O.lpc), supervirial 
(Qvir = T5) Plummer sphere in Fig. 7(b) 



observed distribution (the solid line in Fig. 7(b) I; however, 
the cluster rapidly expands, and after 1 Myr the cluster dis- 
tribution sits on the observed distribution (the dashed line). 
Following further evolution and expansion of the cluster, the 
E distribution lies to the left of the observed distribution af- 
ter 5 Myr (the dashed-dotted line) and 10 Myr (the dotted 
line). 

We note that in both these evolutionary scenarios the 
clusters retain the same Plummer sphere morphologies (see 
e.g. Fig. |l(a)[ ) throughout their dynamical evolution, and 
their virial ratios are essentially constant. 



3.2.2 Fractal clusters 

Recently, it has been postulated that the Orion Nebula Clus- 
ter (ONC) may have been initially subvirial, and under- 
gone 'cool' collapse, leading to dynamical mass segregation 
within 1 Myr and the formation of Trapezium-like s ystems 
(|Allison et al.ll2009bl . l20Tb1 ; lAllison fc Goodwin|[201ll ). Such 
a scenario is also c onsistent with the observed binary separa- 
tion distribution (|Parker et al.|[201ll ). This scenario favours 
the formation of stars in the cluster with a high degree of 
substructure, as a radially smooth cluster cannot mass seg- 
regate with in the age of the ONC, e ven with subvirial initial 
conditions (|Bonnell fc D avics 199^). 

We show the evolution of the E distribution for a highly 
substructured cluster undergoing cool collapse in Fig. [8] The 
cluster has a radius of 1 pc, fractal dimension D — 1.6 and 
a 'cool' virial ratio (Qvir = 0.3). In panels (a) and (b) we 
show the morphology of a typical cluster at (a) OMyr and 
after 5 Myr of dynamical evolution (b). In a cluster under- 
going this dynamical scenario, the stars in the clumpy areas 
of the fractal interact, and the cluster collapses to a central 
concentration after 1 Myr. This central concentration has 
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Figure 8. The typical dynamical evolution of a fractal (D = 1.6) cluster undergoing cool collapse (virial ratio Q v i r = 0.3). In panel (a) 
we show the morphology before evolution, and in panel (b) we show the morphology after 5 Myr. In panel (c) we show the E distribution 
at OMyr, 1 Myr, 5 Myr and 10 Myr by the solid, dashed, dash-dotted and dotted lines, respectively and the data from Bressert et al are 
shown by the red line. In panel (d) we show the E distribution from the simulation in panel ( c) at 5 Myr (the dot-dashed line ) compared 
to the other 9 simulations at 5 Myr. In panel (e) we show the evolution of the Q-paramctcr ( Cartwright & Whitworth 2004) during the 
first 10 Myr of this cluster's evolution (the dot-dashed line, and all other simulations are also shown for comparison). Finally, we show 
the average evolution of the virial ratio in the simulations (with 1-ct uncertainties) in panel (f). 
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a similar appearanc e to a Plummer sphere or King profile 
ison et al.l 120101 ) ■ and the cluster then relaxes and ex- 
pands over the following lOMyr. Note that t he cluster ob- 
tains a maximum central density after ~1 Myr (jParker et al.l 
l201ll ), but Fig. |8(c)| shows that the distribution of E has 
already moved to lower densities compared to the distribu- 
tion at OMyr. This is due to ejections from the substructure 
before the collapse of the cluster, which shifts the E distri- 
bution to lower densities. 

In Fig. 8(c) the initial E distribution is shown by the 
solid line. After 1 Myr the cluster has already reached its 
densest phase and is now expanding, and the distribution 
of E moves to the left (the dashed line). As the cluster dis- 
solves the E distribution widens, and we show the distri- 
butions after 5 Myr (the dot-dashed line) and 10 Myr (the 
dotted line ). The data for all c lusters in the Solar Neigh- 
bourhood (Brc ssert et aljfeoioT ) is shown by the solid red 
line. Fractal clusters evolve very stochastically (|Allison et al.l 
l20ld ; iParker fc Goodwinll2012l ) and there is a variation of 
E due to this. We show the spread in E at 5 Myr in our 
10 simulations in Fig. 8(d) We show the particular clus- 



10 Myr of evolution. The spread in E at 5 Myr between the 



10 simulations is shown in Fig. 9(d) 



ter in Figs. 8(a) 8(c) by the thick dot-dashed line, and the 
remaining simulations by the coloured lines. 

We also show the evolution of the Q-parameter in 
Fig. |8(c)| Again, the cluster in Figs. |8(a)||8(c")| is shown by 
the thick dot-dashed line. We choose this cluster as a 'typi- 
cal' example because neither the E distribution nor the evo- 
lution of the Q-parameter is an outlier, though one could 
also choose other examples which fit this criteria. The ini- 
tially substructured cluster has Q < 0.5, and the dynamical 
evolution in our simulation quickly erases this substructure 
within the first 0.5 Myr before reaching a large (> 1) Q- 
parameter, indicating a centrally concentrated geometry. 

Finally, we show the evolution of the virial ratio for 



this collapsing, substructured cluster in Fig. 8(f) The clus- 
ter begins subvirial (Q v ir = 0.3), but rapidly reaches virial 
equilibrium (Q v ir = 0.5) in the first 1 Myr due to violent 
relaxation of the substructure , and the subseque nt global 
collapse of the cluster (see also I Allison et al.ll2010T ). The be- 
haviour of the virial ratio is qualitatively similar for all 10 
simulations (unlike the evolution of the Q-parameter) and 
we plot the average of all 10 simulations, with 1-a uncer- 



tainties, in Fig. 8(f) 



The evolution of a substructured, supervirial (Q vlr = 
1.5) cluster undergoing expansion is shown in Fig [9] We 
show the morphology of the (typical) cluster before dynam- 
ical evolution in panel (a), and the corresponding morphol- 
ogy after 5 Myr of evolution in panel (b). Again, we choose 
this simulation because it is not an outlier in terms of its 
global density (E) or its structure (the evolution of the 
Q-parameter). The cluster has expanded rapidly, but the 
clumps within the fractal have also evolved so that the clus- 
ter resembles a large association without a central concentra- 
tion. This is confirmed by the evolution of the Q-parameter 



(the dot-dashed line Fig. 9(e) I, which remains below the 



substructure cut-off value of 0.8 throughout the cluster's 
evolution. 

We show the evolution of the E distribution for this 
cluster in Fig. 



There is also a wide spread in the evolution of the Q- 
parameter for these expanding, clumpy clusters. Nine simu- 
lations remain substructured (Q < 0.8) throughout the clus- 
ter's lifetime, but in one simulation the substructure evolves 
to a centrally concentrated profile (the (upper) purple line 
in Fig. 9(e) I. There is also some variation in the virial ratio, 
and we plot the average of all 10 simulations, with l-cr un- 
certainties, in Fig. |9(f)| All the clusters remain supervirial, 
but some more so than others, hence the large uncertainties. 

Finally, we note that for both the cluster undergoing 
cold collapse, and the rapidly expanding cluster, the E dis- 
tribution for the fractal widens. Therefore, if stars were to 
form in a fractal distribution, it would be difficult to rule 
out this formation scenario by analysing the shape of the E 
cumulative distribution. 

Note that we would expect a great deal of dynami- 
cal interaction in the substructure as the supervirial cluster 
expands. Even though it does not reach a centrally con- 
centrated dense phase (like the initially subvirial cluster), 
pockets of sub structure are s ufficie ntly dense enough to pro- 
cess binaries ([Parker et al.l l201ll ). and planetary systems 
jParker fc Quanj|2012l 'l. 



9(c) 



The distribution quickly widens, and 
at 5 Myr sits on the observed distribution (the dot-dashed 
line). In this evolutionary scenario, the median E value de- 
creases by over three orders of magnitude during the first and associations. We demonstrate this in Figs. 5(a) 5(c) 



4 DISCUSSION 

In Section we presented the cumulative distribution for 
stellar surface densities, E, for a range of different single clus- 
ter morphologies. We divide the morphologies into two sub- 
groups; radially smooth (the Plummer spheres and King pro- 
files), and substructured (the fractals and associations). We 
have deliberately chosen the radii of the clusters so that the 
median of the E distribution is overlaid upon the median of 
the en semble observational data presented in lBressert et al.l 
(|2010h . 

From inspection, the radially smooth clusters are able to 
reproduce the shape of the observed cumulative distribution, 
and the King profile, which is more centrally concentrated 
than the Plummer sphere, is an especially good fit. On the 
other hand, the clusters that are initially substructured have 
a narrower cumulative distribution than the observed data. 
We note that varying the parameters of these clusters (over- 
all radius, fractal dimension, association node radius) does 
not widen this distribution. 

We have also demonstrated, with the KS test, that only 
one single cluster morphology, a K i ng pr ofile, is consistent 
with the dataset of Bressert et al] |2010h . This is not sur- 
prising, since the observational dataset includes stars from 
different star forming regions, some of which contain sub- 
structure, and others have a more centrally concentrated, 
smooth radial profile. For example, within the dataset are 
Taurus and Ophi uchus. A study of substru c ture u sing the 
Q-parameter by ICartwright fc Whitworthl |2004l ') showed 
Ophiuchus to have a smooth radial profile (Q = 0.85), 
whereas Taurus is substructured, with a low fractal dimen- 
sion (Q_^045V 

iBressert et all |2010l ) point out that the smooth cumu- 
lative distribution in the data is evidence of hierarchical star 
formation, rather than a bi-modal distribution of clusters 
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Figure 9. The typical dynamical evolution of a fractal (D = 1.6) cluster undergoing warm expansion (virial ratio Q v ; r = 1.5). In panel (a) 
we show the morphology before evolution, and in panel (b) we show the morphology after 5 Myr. In panel (c) we show the E distribution 
at OMyr, lMyr, 5 Myr and 10 Myr by the solid, dashed, dash-dotted and dotted lines, respectively and the data from Bressert et al are 
shown by the red line. In panel (d) we show the E distribution from the simulation in pan el (c) at 5 Myr (dot-dashed line ) compared 
to the other 9 simulations at 5 Myr. In panel (e) we show the evolution of the Q— parameter llCartwright fc Whitwortbl l2004) during the 
first 10 Myr of this cluster's evolution (the dot-dashed line, and all other simulations are also shown for comparison). Finally, we show 
the average evolution of the virial ratio in the simulations (with 1— cr uncertainties) in panel (f). 
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where we place an association in a much more sparse back- 
ground. The histogram displays obvious bimodality, and 
the cumulative distribution shows a prominent kink, which 
would be observed if star formation was heavily bi-modal. 
However, we also caution that if this distinction were to be 
blurred, then a E distribution would not be able to detect a 
more diffuse association against a background of field stars, 



as we demonstrate in Figs.[5(d)H5(f) 



We have also shown the effects of dynamical evolution 
of the clusters on the distribution of E. Firstly, we have mod- 
elled a Plummer sphere with a s imilar half-mass radius to 
that on the ONC today (0.8 pc, [Hillenbra nd fc Hartmannl 
1998). Even for a cluster in virial equilibrium, over a pe- 
riod of 10 Myr the distribution of E evolves to lower densi- 



ties as the cluster expands through interactions. In Fig. 7(a) 



we show the initial distribution (the solid black line), and 
the distributions at 1, 5, and 10 Myr (the dashed, dot- 
dashed and dotted lines, respectively). The distributions at 
and 1 Myr are virtually i dentical, and we note that the dis- 
tribution of E observed bv Bres sert et al. (2010) is for young 
(^lMyr) star forming regions. Therefore, if all star forma- 
tion was to occur in virial equilibrium then we would expect 
the distribution of E to be indicative of the primordial den- 
sity of the star forming region for which it was measured. 

Secondly, we have modelled an initially very dense 
Plummer sphere with a half-mass radius of only 0.1 pc, 
but with an initially supervirial ratio (Q v ir = 1-5), so 
that the cluster immediately expands. It has been argued 
that rapid gas expulsion caused by the first supernovae 
or the collective effect of winds from OB stars in a star 
clust e r can cause the expansion of a cluster (e.g. Tutukoyl 
1978; Bastian & Goodwin! liool ; iGoodwin fc Bastianl 1200a 



Bastian et al.l 120081 . and references therein), altho ugh this 
effect may not be as signific ant as first thought (e.g. lBastianl 
l201ll ; iKruiissen et alj|2012h . We show the effect of this sce- 



nario on the E distribution as a function of time in Fig. 7(b) 



The cluster is initially very dense, but rapidly expands so 
that after 1 Myr the density is similar to that obs erved for 
local star forming regions bv lBressert et afl l|2010T l. 

Recently, a dynamical evolution scenario for the ONC 
has been developed, in which the cluster is initially sub- 
structured, and subvirial. The cluster collapses to a cen- 
trally concentrated sphere which appears very similar to a 
King or Plummer profile. This scenario can allow dynam- 
ical mass se gregation of the highest m ass stars within the 
first lMyr ij Allison et al.ll2009bi, l2010ri, and the form ation 
of Trapezium-like systems ( Allison fc Goodwin! l201ll ). Fur- 
thermore, such initial conditions heavily process a primor- 
dial binary population (which is consistent with the ob- 
served serjaxa^ioji dj^tribution and binary fraction in the 
ONC. | Parker et "aTll201lh and can affect planetary systems 
jParker fc Quanj I2012T ). The distribut ion of E for such a 
cluster is shown in Fig. [8] The cluster rapidly collapses, and 
begins expanding after 1 Myr. The increased tail in the dis- 
tribution at low densities is due to ejections from within the 
substructure. Whilst the change in the E distribution with 
time is not as extreme as in the supervirial Plummer sphere, 
it still demonstrates that if a cluster were to evolve in this 
way, the measured E value at 1 Myr is highly unlikely to be 
primordial or indicative of the amount of previous dynami- 
cal processing. 

For completeness, we have also shown the evolution of 



the E distribution in an initially supervirial fractal (Fig[9]). 
The overall evolution is similar to the supervirial Plummer 
sphere, though not as extreme, due to the much lower initial 
density of the fractal resulting in fewer violent ejections. Un- 
like the subvirial fractal, this cluster retains its substructure 
as it evolves to form an association, and the structural differ- 
ence between these clusters is highlighted by the evolution 



of the Q-parameter in Figs. 8(d) and 9(d) 



The measurement of the E distribution for various dy- 
namical scenarios has shown that, at best, it reflects the 
state of the cluster or star forming region now, and in the 
absence of further information should not be used to make 
inferences on the birth environment of stars or planets. Fig- 
ures [7] [5] and [5] demonstrate that the distribution is likely 
to drastically change in the first 5 Myr or so of the cluster's 
lifetime, unless all stars form in a relatively sparse cluster in 
virial equili brium (recent work b y iGieles. Moeckel fc Clarke! 
|2012l ) and lMoeckel et all (|2012h has also pointed out that 
the E distribution evolves to lower densities through dynam- 
ical interactions). 

However, the E distribution used in tandem with the 
Q-parameter and age estimates may provide information on 
the initial density and level of substr ucture. As substruct ure 
is erased quickly through dynamics ( Allison et al ] l2O10l ). if 
the Q-parameter is high (Q > 0.8) we can postulate that dy- 
namical evolution may have occurred and the present den- 
sity may not be primordial. Conversely, as young, unevolved 
regio ns are likely to be substructured (|Sanchez fc A lfaro 
120091 ). then if the Q-parameter is low (Q < 0.8) the the 
E distribution likely reflects the primordial density of that 
cluster or star forming region. 

The volume-limited census of star for ming regions in 
the lo cal Solar neighbourhood presented by iBressert et alj 
(|2010h suggests that most regions would not be dense enough 
to affect the formation and evolution of binary stars and 
planetary systems (especially if these regions formed in virial 
equilibrium). Unfortunately, no data from the central region 
of the ONC is i nclud ed in the sample of stars from Orion in 
IBressert et~ai1 (|2010h due to crowding of stars in the images. 
The contribution of the Orion region in their sample is 2696 
out of a total of 3857 stars. If the central region of the ONC 
was included, then the number of stars from Orion would in- 
crease by ~1000 and it is unclear how this would change the 
E distribution in the local neighbourhood (though it would 
likely c ontribute to the high-density region of this distri- 
bution, iBate et al.lll998| ). The arguments for a dynamically 
active ONC, presente d bylKroupa e t al. (1999); Scall v et al.l 
(1200511: IParker et all (120091 ); I Allison et all (l2009bl . |2010| ); 
I Allison fc Goodwin! (|201lh : IParker et aLll|201ll ), suggest 
that this cluster underwent dynamical interactions that have 
affected binary stars and planetary systems. 

If one assu mes the uniform clu ster mass function of 
N(M) oc M" u 2 a l|Lada fc Ladal |2003| 1. and also that clus- 
ters are uniformly disrupted (independent of mass), then the 
same mass of stars enters the Galactic field from one very 
high mass cluster, as does stellar mass from many smaller 
clusters with a combined stellar mass equal to the high-mass 
cluster. 

We can postulate that the local Galactic field may 
be populated by the sum of ONC-like, and low-mass star- 
forming regions (if all of these clusters eventually disperse). 
Thus, a different interpretation of the E distribution is 
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that the dominant star forming event in terms of mass 
is the ONC, which contributes 25 per cent of mass to the 
field. If the ONC is dynamically old l|Kroupa et all 1 19991 ; 
lAllison et alll2010l ; IParker et aflboill ) then the stars from 
the ONC with low surface densities today may have been 
in a much more dense region in the past (recall e.g. Fig. [5] 
where a fractal collapses to a central concen tration, ejecting 
a halo of stars in the process |Allisonll2012T )'). 

We therefore suggest that as a lower limit, ~25per 
cent of stars originate in a dense cluster, but primor- 
dial substructure in star-forming regions can potentially 
disrupt binaries and planets without the need for hig h 
global densities (|Parker et al.ll201ll ; [Parker fc Quanzll2012l ). 
This substructure could increase this fraction of dynami- 
cally susceptible systems to ~50per cent. Such a scenario 
would then have im plications for the formation and sur- 
vival of binaries (e .g. lKroupa et al.ll 19991 ; [Marks et aL 20 111: 
Parker et al.l 201 ll) and p l anetary systems feonnell et al. 



2001; Scallv & Clarke 2001; Adams ct al. 2006tToiczak et al 



200sl ; IParker fc Goodwinll2009l ; IParker fc Quanzll2012l ) 

Finally, we suggest that in order to determine both the 
current and initial density, and structure of a star forming 
region, a combination of tools must be used. A low, global 
density does not necessarily preclude the possibility of lo- 
cally dense substructure, which can be identified using the 
Q-parameter. Similarly, a low global density now does not 
rule out the possibility than the cluster was originally sev- 
eral orders of magnitude more dense at an earlier time. In 
the near future, accurate measurements of the positions and 
velocities of stars using data from the GAIA mission and 
its associated spectroscopic surveys (e.g. lRandicbll2012l ). in 
tandem with numerical simulations, will enable a detailed 
dyn amical history of nearby star forming regions to be made 
(e.g. lAllisonll2012l ). 



5 CONCLUSIONS 

In this paper, we have modelled both static and evolving 
star clusters with different morphologies, and determined 
the distribution of stellar surface densities, E, and the Q- 
parameter. In our A^-body simulations we determine the E 
distribution at 0, 1, 5 and lOMyr. Our conclusions are the 
following: 

(i) The distribution of E is degenerate. Many different 
morphologies reproduce a smooth, continuous distribution 
of stellar densities. We find that a substructured star form- 
ing region (either a fractal, or an association) has a nar- 
rower cumulative distribution than a centrally concentrated 
spherical cluster, such as a Plummer sphere or a King pro- 
file. However, these morphologies can be much more easily 
distinguished by other methods , such as the Q-parameter 
jCartwright fc Whitwortnl 120041 1 . Characterizing substruc- 
ture in clusters is important, as this has recently been shown 
to facilitate dynamical mass segregation and disrupt binary 
and planetary systems in young clusters IIAUison et al .1120101 : 
IParker et al.ll201ll : IParker fc Quanzll2012! )~ 

(ii) No single morphology is a good match to th e ob- 
served E distribution presented in lBressert et al.l (|2010h . We 
interpret this as being due to the superposition of E values 
from the different star forming regions in the observational 
sample. 



(iii) The E distribution can be significantly shifted to 
lower densities for both substructured and smooth clusters 
due to early dynamical evolution of the cluster. In the case 
of a supervirial (expanding) cluster, the median E value can 
be similar to the observed distribution after just 1 Myr of 
dynamical evolution. The dynamical evolution of a fractal 
causes the distribution of E to widen, as well as shifting it 
to lower densities. Even an initially dense cluster in virial 
equilibrium will expand, due to two-body relaxation and 
push the E distribution to lower values in the first 10 Myr 
of the cluster's life. 

(iv) Substructure is quickly erased through dynamical 
interactions, and therefore the Q-parameter reaches high 
values rapidly. If all star forming regions start substructured, 
then a low Q-parameter could indicate that the E distribu- 
tion reflects the primordial density of the region. This local 
substructure could facilitate dynamical processing of sys- 
tems, without the need for high global densities. 

(v) Although the dynamical scenarios that drastically 
alter the E distribution (a collapsing fractal cluster, or su- 
pervirial, initially dense Plummer sphere), can be viewed 
as extrema in terms of the formation and evolution of clus- 
ters, they represent possible initial conditions for the ONC 
llKroupa et al.ll 19991; (Parker et alfeoogl ; lAllison et aTll2009bl . 
|2010| ; IParker et al.l 120111 ). which would contribute ~25per 
cent of stars to the local volume-limited E distribution if it 
were to be included in the sample. If primordial substruc- 
ture is taken into account, the fraction of stars that could be 
affected by dynamics in a cluster could be as high as 50 per 
cent. 

(vi) If star forming regions form in virial equilibrium, 
then dynamical relaxation cannot significantly alter the E 
distribution. In such systems, the E distribution is a tracer 
of the initial density of a star forming region. 

The combination of points (iii), (iv) and (v) lead us to 
suggest that the E distribution in local star forming regions 
may indicate a relatively quiescent natal environment for 
star formation, but the possibility of a more dynamically ac- 
tive initial environment for up to 50 per cent of the stars can- 
not be ruled out. We suggest that the E distribution in local 
star forming reg i ons should not be over-interpreted (see also 
iBate et al.ll 19981 ; iGieles et al.ll2012l ). and that one should use 
a range of different metrics (e.g. a minimum spanning tree 
analysis such as the Q-parameter in tandem with age es- 
timates, KS_tests and "inverse binary population synthesis" 
|Kroupalll995al lbhl to determine the past and present state 
of star forming regions. Detailed kinematic information from 
the GAIA mission and its associated spectroscopic surveys 
will enable a detailed dynamical history of field and cluster 
stars to be made, which could in principle be used to infer 
the density of the birth environment of most stars. 



ACKNOWLEDGEMENTS 

We thank the anonymous referee for their comments and 
suggestions, which improved the original manuscript. We 
also thank Eli Bressert and Nate Bastian for kindly provid- 
ing us with the cumulative distribution from iBressert et al.l 
(2010), and for a thorough critique of an earlier draft of this 
work. Thanks also to Richard Allison for useful discussions. 



© 2012 RAS, MNRAS 000,rTHT5l 



Characterizing the state of clusters 15 



The simulations in this work were performed on the BRUTUS 
computing cluster at ETH Zurich. 
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